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Abstract. The Riemannian lOj symbols are spin networks that assign an amplitude 
to each 4-simplex in the Barrett-Crane model of Riemannian quantum gravity. This 
amplitude is a function of the areas of the 10 faces of the 4-simplex, and Barrett and 
Williams have shown that one contribution to its asymptotics comes from the Regge 
action for all non-degenerate 4-simplices with the specified face areas. However, we 
show numerically that the dominant contribution comes from degenerate 4-simplices. As 
a consequence, one can compute the asymptotics of the Riemannian lOj symbols by 
evaluating a 'degenerate spin network', where the rotation group SO (4) is replaced by 
the Euclidean group of isometries of M^. We conjecture formulas for the asymptotics of 
a large class of Riemannian and Lorentzian spin networks in terms of these degenerate 
spin networks, and check these formulas in some special cases. Among other things, this 
conjecture implies that the Lorentzian lOj symbols are asymptotic to 1/16 times the 
Riemannian ones. 



1. Introduction 

In the Ponzano-Regge model of 3-dimensional Riemannian quantum gravity an am- 
plitude is associated with each tetrahedron in a triangulation of spacetime. The amphtude 
depends on the tetrahedron's six edge lengths, which are assumed to be quantized, taking 
values proportional to 2j + 1 where j is a half-integer spin. One can compute this ampli- 
tude either by evaluating an SU(2) spin network shaped like a tetrahedron, or by doing an 
integral. Approximating this integral by the stationary phase method, Ponzano and Regge 
argued that when all six spins are rescaled by the same factor A, the A — s- oo asymptotics of 
the amplitude are given by a simple function of the volume of the tetrahedron and the Regge 
calculus version of its Einstein action. Nobody has yet succeeded in making their argument 
rigorous, but Roberts ^ recently proved their asymptotic formula by a different method. 
This result lays the foundation for a careful study of the relation between the Ponzano-Regge 
model and classical general relativity in 3 dimensions. 

Our concern here is whether a similar result holds for the Barrett-Crane model of ^- 
dimensional Riemannian quantum gravity Q . In this model an amplitude is associated with 
each 4-simplex in a triangulation of spacetime. This amplitude, known as a lOj symbol, 
is a function of the areas of the 10 triangular faces of the 4-simplex. Each triangle area 
is proportional to 2j + 1, where j is a spin labelling the triangle. The amplitude can be 
computed by evaluating an SU(2) x SU(2) spin network whose edges correspond to the 
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triangles of the 4-simplex: 




J45 



There is also an integral formula for the lOj symbol 0]. The problem is to understand the 
asymptotics of the IQj symbol as all 10 spins are rescaled by a factor A and A oo. 

Barrett and Williams |^ applied a stationary phase approximation to the integral for the 
lOj symbol, focusing attention on stationary phase points corresponding to nondegenerate 
4-simplices with the specified face areas. They showed that each such 4-simplex contributes 
to the lOj symbol in a manner that depends on its Regge action. They pointed out the 
existence of contributions from degenerate 4-simplices, but did not analyse them. 

In Section 2 of this paper we begin by applying Barrett and Williams' estimate of the 
lOj symbols to the case where all 10 spins are equal. We find that the contribution of 
their stationary phase points to the 10 j symbols is of order A~^/^. However, our numerical 
calculations show that the 10 j symbols are much larger, of order A^^. This means we must 
look elsewhere to explain the asymptotics of the lOj symbol. 

In Section 3 we analyse the contribution of 'degenerate 4-simplices' to the integral for the 
lOj symbols. They do not correspond to stationary phase points in the integral for the lOj 
symbols; instead, the integrand has a strong maximum at these points. We argue that the 
contribution of a small neighborhood of these points is asymptotically proportional to X~^. 
We give a formula expressing the constant of proportionality as an integral over the space 
of degenerate 4-simplices. We also reduce this to an explicit integral in 5 variables. 

In Section 4, we numerically compare these results to the lOj symbols as calculated using 
the algorithm developed by Christcnsen and Egan Our formula for the contribution of 
degenerate 4-simplices closely matches the actual asymptotics of the lOj symbols. Thus, 
even though our argument that these asymptotics are dominated by degenerate 4-simplices 
is not rigorous, we feel confident that the resulting formula is correct. 

In Section 5 we discuss a new sort of spin network, associated to the representation 
theory of the Euclidean group, which arises naturally in our analysis of the contribution 
of degenerate 4-simplices. Generalizing our results on the Riemannian lOj symbols, we 
conjecture formulas for the asymptotics of a large class of Riemannian spin networks in 
terms of these new 'degenerate spin networks'. We verify this conjecture in a number of 
simple cases. 

In Section 6 we formulate a similar conjecture for Lorentzian spin networks. Taken with 
the previous one this conjecture implies that the A ^ oo asymptotics of a Lorentzian spin 
network in this class are the same, up to a constant, as those of the corresponding Riemannian 
spin network. For example, as A — > oo, the Lorentzian lOj symbol should be asymptotic 
to 1/16 times the corresponding Riemannian 10 j symbol! We conclude by presenting some 
numerical evidence that this is the case. 

2. Stationary Phase Points 

The lOj symbols can be defined using a Riemannian spin network — also known as a 
'balanced' spin network |^ — whose underlying graph is the complete graph on five vertices. 
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The ten edges of the graph are labelled with half-integer spins jki = 0, ^, 1, . . . , where k and 
I refer to the vertices connected by each edge. In this approach, an edge labelled by the spin 
j corresponds to the representation j (g) j of Spin(4) = SU(2) x SU(2), and the IQj symbols 
are computed using the representation theory of this group. 

However, to analyse the asymptotics of the lOj symbol, it is easier to use the integral 
formula due to Barrett j^]. This is: 




■^y'' > k<i 



where is the unit 3-sphere in equipped with its usual Lebesgue measure, (j)ki is the 
angle between the unit vectors hk and hi, and the kernel is given by: 

(2) i^f (</>) := 



smq) 

The normalizing factors in the integral come from the fact that the volume of the unit 
3-sphere is 27r2. 

The lOj symbol gives the amplitude for a 4-simplex with specified triangle areas. Each 
vertex of the above graph corresponds to a tetrahedron in this 4-simplex, and each edge 
of the graph corresponds to the unique triangle shared by two of these tetrahedra. In this 
picture, the spin jki determines the area of the triangle shared by the kth and /th tetrahedra. 
The precise formula for this area is somewhat controversial ]lO| , [Tl|] , but given the integral 
formula for the lOj symbols, we find it convenient to assume the area is proportional to 
2jfc( + 1. Ignoring the constant factor, we thus define triangle areas by: 

a-ki = '2 jki + 1- 

In what follows, we study the behaviour of the lOj symbol as all these triangle areas aki 
are multiplied by a large integer A. This is not the same as multiplying the spins jki by A. 
However, note that as jki ranges over all spins, aki ranges over all positive integers. This 
means that if we multiply the areas Uki by any positive integer A, we can find new spins Jki 
corresponding to the new areas by solving Xajk — 2 Jki + 1. 

It is shown in |^ that the integral in (^) is nonnegative. Accordingly, we will concentrate 
on analysing the asymptotics of the absolute value of the lOj symbol, given by the integral 
alone: 



(3) 




'(S^)' k<l 

where for simplicity we have left out the spins labelling the spin network edges, which are 
now Jki- 

Barrett and Williams |^ express the numerator in the kernel as a difference of expo- 
nentials, which allows them to rewrite the integrand in equation (|^) as a sum of 2^*^ terms, 
each consisting of an exponential times the function: 

(4) /(/ii,...,/i5)=n-V- 
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This function is unbounded as any of the (jjki tend to zero or vr, but if the regions of the 
domain where that occurs are set aside for a separate analysis, the integral over the remaining 
region becomes amenable to a stationary phase approximation [Q. The relevant phase is 
simply: 

(5) S{hi, . . . , fag) ^ Xaki(t>ki- 

k<l 

The phases of the 2-*^" exponentials arising from the product of the kernels include variants of 
S with all possible signs for the ten terms. However, Barrett and Williams use the invariance 
of the integrand under the transformations hk — > — /ifc to break the integral into 2^ identical 
parts, according to whether hk — ±nfc, where nk are the outward normals of a 4-simplex 
whose five tetrahedra lie in the hyperplanes normal to the hk- They then focus on the 
region where all the hk = rik, and show that of the 2-'^° phases, only S and its opposite 
have stationary points in this region, and they occur when the cj)ki are the angles between 
the outward normals to the tetrahedra of a 4-simplex whose ten faces have areas given by 
Xaki- In this case the hk can be interpreted as the outward normals to the five tetrahedra, 
the (pki are the angles between these normals, and Xau is the area of the face shared by the 
tetrahedra numbered k and I. Interpreted this way, S is precisely the Regge action for the 
4-simplex. 

This means that the relevant stationary phase points are those of a much simpler integral: 



(6) ^ / fihu . . . , /^5)(e'^^'--''=) + e--(''--''=)) §i • • • 

To carry out a stationary phase approximation of this integral we must note that the sta- 
tionary phase 'point' determined by the geometry of a 4-simplex is not actually a single 
point in the 15-dimensional manifold (S^)^, but rather a whole 6-dimensional submanifold, 
since the geometry of the 4-simplex is invariant under the 6-parameter rotation group SO (4). 
However, the same invariance can be used to remove this complication. Since the functions / 
and S are invariant under the action of S0(4), they pass to the quotient space (S^)^ / S0(4), 
and we obtain: 

(7) :^7;^ / fix){e^'^^^ + e-^^(-)) 
2(2tt)^" 7(53)5/ so(4) 

where d/i is the result of pushing forward Lebesgue measure on (S^)^ to this quotient space. 
A factor of 1/2 here accounts for the complication that simultaneously negating all the hk 
has the same action as — /, an element of SO (4), and so the discrete symmetries used to 
obtain (||) already entailed taking the quotient by a 2-element subgroup of SO (4). 
Now, the stationary phase approximation [ p^ of an n-dimensional integral 

/(x)e'^(^) dx 

is a sum over stationary points Xi of the function S: 

(2-)"/^y:^-J^exp 



ZTT 

iSix,) + —a{H{x^)) 



where H{xi) is the matrix of second partial derivatives of S at the point a:^, and a{H{xi)) 
is the signature of this matrix, i.e., the number of positive eigenvalues minus the number of 
negative eigenvalues. This approximation assumes there are finitely many stationary points, 
all with det H{xi) ^ 0. Applying this to the case at hand, and neglecting points where some 
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of the angles are or tt, we obtain Barrett and Williams' stationary phase approximation 
of the 10 j symbols: 



(8) 




1 



E 



(27r)ii/2 |detiJ(a;,)|i/2 



S{xi 



where the sum is over the 4-simplices Xi with the specified face areas. 

We can easily say something about the A — *■ oo behavior of this quantity without actually 
evaluating it. Each stationary point Xi is independent of A, so f{xi) will be constant as 
A — !■ oo, while S{xi) will grow linearly, as will its matrix of second derivatives, H{xi). This 
is a 9 X 9 matrix, since {S'^)^ / S0(4) is 9-dimensional, so its determinant will grow as A^. It 
follows that: 



(9) 




0(A 



-9/2N 



However, cancellation between different stationary points could in principle make this a mis- 
leading overestimate. Thus it seems worthwhile to explicitly evaluate the sum in equation (H), 
at least in a simple special case. 

In what follows we evaluate this sum for the case of a lOj symbol with all edges labelled 
by the same spins. In other words, we calculate the stationary phase approximation of the 
amplitude for a regular 4-simplex. 

The first step, which would be useful more generally, is to find explicit coordinates on 
the quotient space {S'^)^ / S0(4) and describe the measure dji in these terms. To do this, we 
exploit the fact that any configuration of the five unit vectors can be rotated into one 
that belongs to a 9-dimensional subspace of the original domain. Specifically, if we express 
the vectors hk in polar coordinates 

hk = (cosV^fc, sinV'fc cos6'fc, sin V'fc sin 6*^- cos , sin V'fc sin 6*^- sin 0^ ) , 

where < "^fci ^fc < and < 0^ < 27r, any set of hk can be rotated in such a way that the 
following restrictions are met: 

(10a) /ii = (1, 0, 0, 0) Vi=fi = '/'i=0 

(10b) h2 = (cosV-a, sinV'2, 0, 0) 82 = (t>2 = 

(10c) /13 = (cost's, sin cos6'3, sink's sin^s, 0) 03 = 0. 

By means of this 'gauge-fixing' we can use the remaining 9 variables "03: V'4j "05 j 
^4, ^5, (/)4, (j)^ as coordinates on the quotient space. To describe the measure in these 
coordinates recall that in polar coordinates, Lebesgue measure on the unit 3-sphere is given 
by 

(11a) dhk = sin^Vfe sinflfe dV'fc dO^ d(l)k- 

Since hi is completely fixed we omit dhi from the formula for d/i, only inserting a factor 27r^ 
due to the volume of the 3-sphere. Since /i2 and /13 are partially fixed we replace dh2 and 
dh^ by 



(lib) 
(11c) 

We thus obtain 
(12) 



dh2 ~ 47r sin ip2 dip2 

dhz = 27rsin^'(/'3 sin03 dip^ dOz 



dfj, — 2n dh2dhzdhidh^. 
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Next, we determine the stationary points of the function S in the case where aU ten spins 
are equal, neglecting points where some of the angles (jjki are or tt — that is, where some 
of the vectors hk are parallel or anti-parallel. Recall that these stationary points correspond 
to nondegenerate 4-simplices having all 10 face areas equal. One obvious candidate is the 
regular 4-simplex. However, we must rule out the possibility that there are other, non-regular 
4-simplices for which all the faces have identical areas. Since we are excluding degenerate 4- 
simplices from the current analysis, we can appeal to a theorem of Bang jist which states that 
if all the faces of a non-degenerate tetrahedron have the same area, they are all congruent. It 
follows that if all ten triangles in a non-degenerate 4-simplex have the same area, they too are 
all congruent. Now, each of the ten edges of a 4-simplex is shared by three triangles, so we 
can treat each edge as a triple of congruent line segments that happen to be superimposed, 
giving a total of 30 in all. If all the triangles were congruent isoceles triangles, each with one 
side of length L, then ten of these 30 line segments would be of length L. However, there is 
no way to partition ten line segments into congruent triples. The same argument rules out 
scalene triangles. So the faces must be equilateral, and the 4-simplex must be regular. 

We must therefore find all sets of unit vectors hk which satisfy the gauge-fixing conditions 
in equation ( |l0|) and form the outward normals of a regular 4-simplex. One choice is: 



m = (1,0,0,0) 

n2 = (-i4^,0,0) 

«3 = (-i-4^,V^,0) 




These vectors have mutual dot products of —j. They represent a single point in our 9- 
dimensional domain: 



(13a) ^/;2 = Va = V'4 = V'5 = cos"i(-i) 

(13b) 6I3 = 6I4 = 6*5 = cos-i(-i) 

(13c) 04 = _cos-i(-5) = 

(13d) 05 = cos-i(-l) = f ■ 



The only other choice comes from interchanging 04 and 05. This yields another regular 
4-simplex. 

Now we take the phase (^ and specialise to the case where all the spins jki are equal, say 
to j. Setting a — 2j + I, and working in our chosen coordinates, we obtain 

(14) 5(V'2, l/'3, '04, V'5, 6*3,^*4, 615, 04, 05) = Aa^COS"^(/lfe • hi) 

k<l 

The partial derivatives of S are all zero at the point described by ([l^). Symbolic computer 
calculations show that at this point the matrix of second derivatives of S has 5 positive 
eigenvalues and 4 negative ones, and a determinant of: 

y^.|(|)"(Aa)^. 

Applying this result to (^ at the point described by (|l3|), working in our chosen coordinates, 
redefining the function / to include the measure (^2|) as well as the kernel denominators (||) , 
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and multiplying by a factor of two to account for the second point where 4>4, and 05 are 
interchanged, we obtain: 



(15) 




57r2 



cos 



10cos"V-i) Aa+ T 



The above expression consists of an oscillating term times a function proportional to A 
In short, the estimate 




at 

is sharp, at least for the case of 10 equal spins. 

However, our numerical calculations of the lOj symbols exhibit very different behaviour. 
Rather than being of order A~^/^, they appear to be much larger, of order A~^. Also, 
they exhibit no discernable oscillations. In the next section we argue that these results are 
explained by the contribution of 'degenerate 4-simplices'. 



3. Degenerate Points 



The absolute value of the kernel 
(16) 



sm aq 
sinc^ 



has maxima when the angle (j) equals or tt, and these maxima become ever more sharply 
peaked as a — > oo. This means that as A ^ oo, the integrand in equation (^ becomes very 
large at 'degenerate points', where some of the vectors hk are either parallel or anti-parallel. 
We have made a detailed study of the 'fully degenerate' points, where all the vectors hk are 
either parallel or anti-parallel. It is plausible to expect the contribution from a neighborhood 
of these points to dominate the integral for the lOj symbol, at least asymptotically as A — > oo, 
since in this region we are integrating a product of kernels, all of which are near their greatest 
possible absolute value. In the next section we shall present numerical evidence that this is 
in fact the case. 

The value of the kernel is positive at (f) — 0, but at (j) — tt its sign is positive when 2j 
is even and negative otherwise, where a = 2j + 1. At first glance, this appears to allow 
the possibility that parallel and anti-parallel degenerate points might cancel each other for 
certain values of the spins. In fact, cancellation can occur only when the lOj symbol vanishes. 
The integrand in ([^) is the product of ten kernels, so it will be positive when all the hk are 
parallel. If hi is replaced by its opposite, the four 0i; will change from zero to tt, leading to 
an overall sign change of (— I)20'i2+ii3+ii4+ji5)_ However, for the lOj symbol to be non-zero, 
the spins at each vertex must sum to an integer [Q. Therefore, if the lOj symbol is non-zero, 
the sign of the integrand will remain positive. The same argument applies if any subset of 
the hk are reversed. 

In the integral in equation (^ , fully degenerate points occur in 3-dimensional submanifolds 
of the domain, but we can apply the same gauge-fixing principles as we used to analyse the 
stationary points, in this case simply fixing hi = (1, 0, 0, 0). There are then 16 discrete 
fully degenerate points in the restricted domain: those where hk = ±hi for k — 2, . . . , 5. 
Restricting the integral to the vicinity of these points we obtain 



(17) 




:= 16 



R , 



dh2 



27r2 



dh^ 
27r2' 



dcg " k<l 

where U is a small open ball around the point {hi, hi, hi, hi) G (S^)'^. 
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We approximate this quantity by noting that when (j) is small, the kernel in ( |l^ ) is close 



to: 

(18) 



sm ac 



We call this quantity the 'degenerate kernel'. As we shall see in Section this is the analog 
of the original kernel in a spin-network formalism where the space of constant curvature, S^, 
is replaced by three-dimensional Euclidean space, and the group SO (4), the isometry group 
of S^, is replaced by the Euclidean group of isometrics of M.^. 

Similarly, the integral over a small subset of (S^)"^ can be approximated by an integral 
over a subset of (K'^)'*, and the angle (pki between unit vectors hk and hi in replaced by 
the Euclidean distance rki = |a:fc — a;;| between vectors Xk and xi in R^. In terms of these 
new Euclidean variables, the restriction hi = (1, 0, 0, 0) is replaced by xi = (0, 0, 0). 

Thus, when A is large, we have: 



(19) 




16 / ll^LM 

•'^ k<l 



dX2 
^l\) TT^ 



2^' 



where now J7 is a small open ball around the origin of 

The approximation ( pj[ ) exhibits very simple scaling behaviour, 
degenerate kernel (|l^) obeys the identity: 



First, note that the 



(20) 



This allows the scaling of (U 



r 



= AK,^(Ar) 

to be deduced from a linear change of variables, 



Vk 



Xxk 



16 



(21) 



16A- 



xu 



k<l 



dxo 



( 



Wk - yi\ 



k<l 



27r2 
27r2 



da; 5 
'2^ 

dys 
27r2' 



where XU is the result of rescaling the neighborhood [/ by a factor of A. If the integral on 
the right hand side of this equation converges in the limit A — > oo, wc obtain the asymptotic 
formula: 



(22) 




16A- 



dog 



k<l 



Ifc! ' 



Wk 



IN 



dy5 

27r2' 



We call this integral the 'degenerate lOj symbol'. Apart from the factor of IGA^^, the integral 
here can be interpreted as the evaluation of a spin network with edges labelled by unitary 
irreducible representations of the Euclidean group: the group of isometrics of Euclidean 
3-space. We discuss this interpretation in more detail in Section ^ 

While elegant, the integral in (|2^) is difficult to compute numerically, very much like the 
integral for the Lorentzian lOj symbol. To obtain a more convenient form for evaluation, we 
make use of the Kirillov trace formula: 

sinla;! 1 



47r 



exp(za; • ^) d£_, 



where a; is a vector in , and the integral is over the unit 2-sphere with its standard measure. 
This formula is easily confirmed by choosing spherical coordinates such that the z-axis is 
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parallel to the vector x. For our purposes we shall rewrite it as follows: 
(23) Ky\x\)= I exp(*x.C) dt 

JS{a) 

where S{a) is the 2-sphere of radius a embedded in , but where is the induced Lebesgue 
measure divided by 47ra, to hide some annoying constants that would otherwise appear in 
this formula. Using this we can rewrite ( p2[ ) as: 

(24) 

16A-2 




dog 



(27r2)4 

7^8 



4 Jx 

2ol2„4 / rl2 



exp(i y^ivk - yi) ■ £,ki) d£,i2 ■ ■ ■ d^45 ^2/2 • • • dy^ 

exp(i(y2, 2/3, y4, J/s) • F{ii2-, ■ ■ ■ £,45)) d^u ■ ■ ■ d^45 dy2--- dy^ 
,£45)) d£i2 ■ --dUs, 



Ix 



where 



k<l 



is a Cartesian product of 2-spheres with the measures described above, and F: A" — > R-'^^ is 
defined by: 

F{£,12, ■••,■^45) = (^^12+^23+ ?24 + £.25, 

— £13 — £,23 + £34 + £35 , 

— £14 — "^24 — £34 + £45 , 

— £15 — £25 — £35 — £45)- 

If 

M^{£€X: Fi£) = 0} 
then the final integral in ( p^ ) will be well-behaved so long as at each point ^ e TV the 
differential of F has maximal rank, namely 12. When this is the case Af is an 8-dimensional 
submanifold of X, and the integral reduces to: 



(25) 




\Jm-' d£ 



Here d£ is the Lebesgue measure on Af induced by the Riemannian metric on X , but divided 
by a factor of Y[k<i ^^rafci, since we have divided the Lebesgue measure on each sphere by 



a factor of 'iiraki- If we choose local coordinates (x^, 



^20 



) on X near £ £ M such that 



x^, . . . , a;^ are zero on A/", then \J{£)\ is the Jacobian determinant of F as a function of the 
remaining 12 variables a;®, . . . , x^^ . 

One way to get solutions of F = is to start with 4-simplices in R'^: that is, 5-tuples of 
points in R'^, together with the ten triangles and five tetrahedra determined by these points. 
Given such a 4-simplex, let £ki be the vector that is normal to the triangle shared by the 
A:th and /th tetrahedra, and has length equal to the area of this triangle. As shown in 
the four vectors £ki normal to the triangles in any one tetrahedron must sum to zero, with 
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appropriate signs: 

-62+63+64+65 = 

—63 — 63 + 64 + 65 = 
(26) -64-64-64+65 = 

-65 - 65 - 65 - 65 = 
62 + 63 + 64 + 65 = 0. 
Each vector appears twice in these formulas, with opposite signs, since the outwards-pointing 
normal to one tetrahedron is the inwards-pointing normal to another. The first four equations 
say that F = 0; the last is an algebraic consequence of the rest. If |6i| = o,ki, the vectors 6i 
thus determine a point in M. 

This suggests that we think of points of J\f as 'degenerate 4-simplices'. However, not 
every point of Af comes from a 5-tuplc of points in M"^ this way. To see this, note that two 
5-tuples in E'^ determine the same point in Af if they are translates of each other. The space 
of 5-tuples modulo translation has dimension 3 x 5 — 3 = 12, but the space of solutions of 
equation ( p6| ) has dimension 30 — 4 x 3 = 18. Thus there are simply not enough 5-tuples 
of points in to account for all solutions of equation (p6|). We shall still call points of Af 
'degenerate 4-simplices', because it is a useful heuristic. However, it is important to take 
this phrase with a grain of salt. 

Of course, for the space Af to be nonempty, there must exist a solution of -F = 0. This 
imposes certain restrictions on the numbers aki- For example, there will be a solution of 

62 + 63 + 64 + 65 = if and only if these 'tetrahedron inequalities' hold: 

ai2 < ai3 + ai4 + ai5 

ai3 < ai2 + ai4 + ais 

ai4 < ai2 + ai3 + ai5 

ai5 < ai2 + ai3 + ai4- 

In general Af will be empty if the four numbers aki corresponding to the faces of any one 
tetrahedron violate the tetrahedron inequalities. In this case the degenerate lOj symbols 
vanish. Similarly, the Riemannian IQj symbols vanish if the spins jki violate the tetrahedron 
inequalities Q , and for the same sort of geometrical reason [|l^ . 

For computational purposes it is helpful to rewrite the integral for the degenerate lOj 
symbols in yet another way. To do this, first we note that we can exploit S0(3) invariance 
to convert ( [25| ) to an integral over the 5-dimensional manifold Af / SO (3), by rotating the £,ki 
in X so that they lie in the 17-dimensional subspace where 63 is fixed at (0, 0, 023), and 64 
lies in the a; > half of the xz-plane. This corresponds to integrating over all possible values 
for 631 inserting a factor equal to the volume of S'(a23) with our chosen measure, which is 
023 ; and also integrating over all possible azimuthal coordinates for ^34 and inserting a factor 
of 2tt. 

In what follows, we will choose coordinates so that the 12 x 12 matrix J(6 whose deter- 
minant we require is block diagonal, with two 3x3 blocks and one 6x6 block. 

It turns out to be convenient to parameterise 64 1 not by its angle from the z-axis, but by 
the length 

si I63 - 64!- 

Since 63+63 = 64+651 these four vectors can be positioned to form the sides of a (possibly 
non-planar) quadrilateral. Then si is the length of one of the diagonals. With the vectors 

63 and 64 fixed and the lengths of 63 ^ind 65 specified, the only remaining freedom this 
quadrilateral has, if it is to remain closed, is the 'hinge angle', ai, between the two triangles 
that meet along the diagonal. The vectors 63 ^^nd 65 have 4 degrees of freedom in all, and 
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specifying ai removes one of them, leaving 3 which break the quadrilateral. With our chosen 
measure on Af, the product of the measure for the coordinates we are integrating over, and 
the Jacobian determinant for the 3 that break the quadrilateral, is: 

1 

(47r)3a23 ■ 

We can treat a second 4-tuple of vectors more or less identically. If the vector ^12 has an 
azimuthal angle of and we define 

S2 ■■= I63 - C12I, 

then the quadrilateral formed by ■^23,^12,^25 and ^24 can be assigned a 'hinge angle' of 
a2- This specifies all the degrees of freedom that allow this quadrilateral to remain closed. 
Once again, the product of the measure for the coordinates we are integrating over, and the 
Jacobian determinant for the 3 that break the quadrilateral, is: 

1 

(47r)3a23 ' 

So far, we have specified 5 degrees of freedom for A/'/S0(3): Si, S2, ai, 0:2 and cj). No 
further continuous degrees of freedom remain. The three vectors we have yet to specify 
must form triangles that complete two quadrilaterals with vectors that have already been 
parameterised, and because these two triangles have a vector in common, there is no 'hinge' 
freedom left. 

Specifically, the three vectors ^14,^15 and ^45 must form a tetrahedron by fitting over a 
triangular base that has, as two of its sides, the vectors: 

V ;= 65 + 65 
w := -^24 - 64- 

Given their fixed lengths, this determines ^14, ^15 and ^45 completely, apart from the freedom 
to locate the apex of the tetrahedron on either side of the plane spanned by v and w. This 
freedom can be accounted for with a factor of 2 in the integral. The final contribution to 
the Jacobian comes from a 6 x 6 block involving all the coordinates of ^14,^15 and ^45, and 
with our chosen measure this is: 

1 

(47r)3 6V{ai4, ais, 045, \w\, \v - w|, l^l) ' 
where V is the volume of the tetrahedron as a function of its edge lengths. 
Combining these results, we can rewrite ( p5| ) as: 

A^^ f f f f f dsi ds2 dai da2 d(p 

967r4a23 Jsi J S2 J 0,2 h V'(ai4, ais, 045, \w\, \v - w\, \v\y 

The integrals over the Si are taken over intervals determined by the four sides of the quadri- 
laterals for which they are the diagonal lengths, and the angular variables range from to 
27r, with the proviso that any part of the domain where the tetrahedron is not geometri- 
cally possible must be excluded. In numerical calculations, this can be dealt with by setting 
the integrand to zero wherever Cayley's determinant formula for the squared volume of the 
tetrahedron yields a negative value. 

We note that the integrand here is unbounded, and we have not proved that ( |27| ) converges, 
but our numerical calculations suggest that it does. 
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4. Numerical Data 

To test our hypothesis that the asymptotics of lOj symbols are dominated by the contri- 
bution of degenerate 4-simphces, we used the algorithm described in to calculate values 
for several sets of lOj symbols. The figure below shows log-log plots for the absolute values 
of the Riemannian lOj symbols as a function of A, where A is the parameter by which the 
areas au were multiplied. The legend shows the base spins jki] multiplying aki by A was 
achieved by replacing the jki with: 

A- 1 

Jki = \]ki H — • 





1/2 
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The lines on the plot show the A ^ asymptotic behaviour described by equation (p7|); the 
integrals were evaluated numerically with Lepage's VEGAS algorithm and found to 
have values of 0.680, 0.341, 0.209, 0.110, 0.0841 and 0.0446 respectively. 

In summary, our numerical data supports the following conjecture. Let us say that the 
spins jki are 'admissible' if for each vertex in the lOj symbol, the spins labelling the four 
incident edges satisfy the tetrahedron inequalities and sum to an integer. Then: 

Conjecture 1. If the ten spins jki are admissible, the A — > oo asymptotics of the Riemannian 
lOj symbols are given by: 




In the next two sections we generalize this conjecture to a large class of Riemannian and 
Lorentzian spin networks, including the Lorentzian lOj symbols. 

The reader may have wondered why we consider asymptotics of the lOj symbols as areas 
are rescaled, instead of spins. The reason is that they are much simpler. We also calculated 
values for sets of lOj symbols where the spins jki were multiplied by A. The figure below 
shows log-log plots for the absolute values of the Riemannian 10 j symbols with spins Xjki, 
where the jki match those shown in the legend for the previous figure: 




Since the areas 2Xjki + 1 are asymptotic to a series of values 2Xjki that are proportional 
to A, the data might be expected to exhibit A~^ scaling. And in fact, in cases A-D, that is 
exactly what is seen. Furthermore, numerically evaluating the integral ( |27| ) with Oki — 2jki 
provided the correct coefRcients for lines on the plot which are asymptotic to the data; these 
coefficients were 2.73, 0.987, 0.472 and 0.165 respectively. 

For cases E and F, the lOj symbols have one or more vertices at the 'border of admissibil- 
ity': that is, vertices where three of the spins labelling incident edges sum to equal the fourth 
spin, making one of the tetrahedron inequalities a strict equality. These vertices are marked 
with heavy dots on the legend. Setting Uki — 2jki in ( p7| ) in these cases yields a coefficient of 
zero, because the domain for at least one variable of integration, the quadrilateral diagonal 
length si, was reduced to a single point. Empirically, the data here appears to scale as A^'^. 
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While lOj symbols with one or two vertices on the border of admissibility have 
asymptotics with spin rescaling, lOj symbols with three or more vertices on the border of 
admissibility appear to decay exponentially as their spins are multiplied by A, as illustrated 
in the log-linear plot below: 




To prove convergence of the partition function in our newly formulated version of the 
Barrett-Crane model Q , it will probably be necessary to understand the A^"^ and exponential 
decay of borderline-admissible lOj symbols under spin rescaling. However, these are more 
delicate phenomena than we are prepared to tackle here. 

5. Degenerate Spin Networks 

Though the Riemannian 10 j symbols are given by an integral over a product of copies 
of S^, we have seen that when we rescale the areas by a large constant A, this integral is 
dominated by the contribution of a very small patch of this space, which can be approximated 
by a product of copies of R^^. Indeed, by a change of variables we can think of the A — > cx) 



ASYMPTOTICS OF lOj SYMBOLS 



15 



limit as one in which the radius of approaches infinity, so that it degenerates to Euchdean 
3-space. 

The Riemannian lOj symbols can also be described using the representation theory of 
S0(4), the isometry group of S^. This suggests that the A — > cx) asymptotics of the lOj 
symbols can be described in terms of the representation theory of the Euclidean group E(3), 
the isometry group of R'^. 

To do this, we introduce certain spin networks associated to E(3) which we call 'degenerate 
spin networks'. These are more closely analogous to the Lorentzian spin networks defined 
in than to the Riemannian spin networks we have been discussing so far, because the 
edge labels are not restricted to discrete values, and the group representations are infinite- 
dimensional. However, all three sorts of spin network form part of a unified theory, as 
outlined below: 



geometry 


signature 


symmetry group 


homogeneous space 


Riemannian 


(++++) 


S0(4) 




degenerate 


(0+++) 


E(3) 




Lorentzian 


(-+++) 


S0(3,l) 





The rotation group S0(4) and the Lorentz group S0(3, 1) consist of linear transformations of 
M"* with determinant 1 that preserve metrics of signature (++++) and ( — respectively. 
Similarly, the Euclidean group E(3) is isomorphic to the group of linear transformations of 

with determinant 1 that preserve the singly degenerate metric with signature (OH — 
The groups S0(4) and S0(3, 1) both 'contract' to E(3), meaning that they have it as a 
limiting case if one forms the isometry group of the metric diag(e, 1, 1, 1) and lets e J, and 
e t 0, respectively. 

This makes it plausible that the asymptotics of not only Riemannian but also Lorentzian 
spin networks can be calculated using degenerate spin networks. In this section we state a 
precise conjecture along these lines for a large class of Riemannian spin networks, and present 
some supporting evidence. In Section^ we do the same for Lorentzian spin networks. 

We begin by describing degenerate spin networks and how to evaluate them. The repre- 
sentations j 1^ j labelling edges of a Riemannian spin network are representations not just of 
Spin(4), but actually of SO (4). As emphasized by Freidel and Krasnov these represen- 
tations are precisely the eigenspaces of the Laplacian on , which is the homogeneous space 
SO(4)/SO(3). Similarly, the representations labelling edges of a Lorentzian spin network 
are the eigenspaces of the Laplacian on hyperbolic 3-space, = S0o(3, 1)/ S0(3) — here 
we must take the connected component of the Lorentz group to get just one sheet of the 
hyperboloid. Following this pattern, the representations labelling edges of a degenerate spin 
network should be the eigenspaces of the Laplacian on = E(3)/ S0(3). 

In fact, there is one representation of this sort for each positive real number a. Any 
complex function on M.^ with V^/ = —a^f can be written as: 

(28) fix) = [ m exp(zC • x) d^, 

where S{a) is the 2-sphere of radius a centered at the origin of M^, and d^ is the induced 
Lebesgue measure divided by Aira. Defining an inner product on these solutions by 

if, 9)= f W)mdi, 
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we form a Hilbert space Ha consisting of all solutions / with (/, /) < oo. This Hilbert 
space becomes a representation of the Euclidean group where each group element g acts via 
{gf){x) — f(g~^x). In fact, this representation is unitary and irreducible. The functions 
exp(i^ • x) form a 'basis', in the sense that any element of TLa can be expressed as in equation 
( psf ) for a unique square-integrable function / on the sphere. 

We define a 'degenerate spin network' to be a directed graph with each edge e labelled 
by a representation of this form, or equivalently, a positive number Og. An intertwiner 
between tensor products of these representations can be defined at each vertex by taking the 
product of functions from the representations labelling the incoming edges, multiplying it by 
the product of complex conjugates of functions from representations labelling the outgoing 
edges, and integrating the result over . Given this, the standard way to evaluate such a 
spin network [|l8| would be to take a 'trace': that is, integrate over a basis label S S{ae) 
for each edge of the graph. The result would be: 



n 



exp(i(a:s(e) - Xt(^e)) ■ Ce) dS,^ , , „ 



ndXy 
0^' 



Here E denotes the set of edges of the graph, V denotes its set of vertices, and the vertices 
s(e) and t{e) are the source and target of the edge e. 

However, just as in the Lorentzian case |l2|, this gives a divergent integral, because the 
integrand is invariant when we simultaneously translate all the vectors Xy e by the same 
amount. More generally, if the underlying graph of our spin network consists of several 
connected components, and we translate the vectors x^ where v lies in any one component, 
the integrand does not change. To keep things simple, let us consider only spin networks 
whose underlying graph is connected. In this case we can sometimes obtain a well-defined 
integral by 'gauge-fixing' one of the vectors rather than integrating over it: that is, setting 
Xy^ = x € M?" for some vertex vi V . If we let V = V — {fi}, this gives the following 
formula for evaluating a degenerate spin network with edges labelled by the numbers Of.: 



(29) I^{a) = / n 



veV 



exp(z(Xs(e) - Xt{e)) ' ^e) d^e 



ndXy 



As in the Lorentzian case one can show that if this integral converges, the result does 
not depend on our choice of the special vertex vi or the point a; G M^. Assuming the integral 
does converge, we can use the Kirillov trace formula (|23| ) to reexpress it as: 



(30) 



When the evaluation of a degenerate spin network converges, it always obeys a very simple 
scaling law as we multiply all the edge labels ae by the same constant A. Using the scaling 
property of the degenerate kernel noted in equation (^), we find that 

(31) /^(Aa) = Al^l-3(IV'l-i) 

where \E\ is the number of edges in the underlying graph and \V\ is the number of vertices. 

We have already argued that the asymptotics of the Riemannian lOj symbols are governed 
by the corresponding degenerate spin network. We can generalize this argument as follows. 
Fix a connected graph. If we label each edge e by a positive integer — or equivalently 
a spin je with Oe = 2je + 1 — we obtain a Riemannian spin network, whose evaluation we 
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define by: 

(32) i^{a)= [ n<(^(^^W'^*(e))) n 

I r\ -TL Z7T 

11 * eeE vev 

Here d{x, y) is the distance between points x, y in the unit 3-sphere in M"* as measured by the 
induced Riemannian metric. This formula is equivalent to the standard integral formula [ pO| , 
except that we have omitted the usual signs in order to simplify the relationship to degenerate 
spin networks. Fixing a small open ball U around some point {x, . . . ,x) S IltjGV" define 
the 'degenerate contribution' to this integral to be: 

(33) /^(a) = 21^1-1 / n^i!(rf(-.(e),-*(e))) n 

eeE veV'^ 



Just as we included a factor of 16 in equation (17), here we include a factor of 21^'"-'^ 
to take into account the contribution of anti-parallel degenerate points; as before there is 
no cancellation between these if the spins labelling edges incident to each vertex sum to 
an integer, as they must for the spin network to have a nonzero value. Using the same 
nonrigorous argument as in Section |^, we see that as A — > oo, 

/«JAa) ^ 21^1-1 / X{KZS\^,^,)-x,^e)\) 

•'^ eeE vev 

(34) .2l-MAl-l-3(m-)/ n<(l^«(e)-2/*(e)l) nil 

•'^•UeeE vev 

= 2l^l-i/^(Aa), 

where U now denotes an open ball around the origin of Hwev ^ ^'^'^ made the change 
of variables ye = \xe- 

In short, this argument suggests that the asymptotics of the degenerate contribution 
to the value of a Riemannian spin network are proportional to those of the corresponding 
degenerate spin network: 

/^(Aa)^ 21^1-1 /^(Aa), 
and we know the latter are very simple: 

/^(Aa) = Al^l-3(l^|-i) jD^^y 
This is particularly interesting when we also have 

/«(Aa)^/<«g(Aa), 

because then we can compute the asymptotics of a Riemannian spin network by evaluating 
a degenerate spin network: 

/^■(Aa)^2l^l-i Al^l-3(l^l-i) /^(a). 

When can we expect this to occur? Clearly we should at least demand that the degenerate 
contribution outweigh the contribution of stationary phase points. A simple power-counting 
argument as in Section^ suggests that the contribution of stationary phase points is of order: 

H^I+3) \V\>2 

(35) /«,t(Aa) = <; 0(A-^) \V\^2 

\V\ = h 
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where the graphs with one or two vertices are different because there is less need for 'gauge- 
fixing'. Comparing these asymptotics to those of the degenerate contribution, we can for- 
mulate the following: 

Conjecture 2. Given a connected graph with more than two vertices and \E\ > ^\V\, or 
two vertices and \E\ > 2, or one vertex and \E\ > 1, as X ^ oo we have 

as long as the integral defining (a) converges and the spins je labelling the edges incident 
to each vertex are admissible. 

Here we say the spins labelling the edges incident to some vertex are 'admissible' if they 
sum to an integer and each is less than or equal to the sum of the rest. We do not yet have 
general criteria for when the integrals associated to Euclidean spin networks converge, and 
as we shall see, the relevant theorems are bound to be a bit different than in the Lorcntzian 
case [ p^ . 

The simplest test of this conjecture is the 'theta network', with two vertices joined by 
three edges, labelled by positive integers a, 6, and c. When the corresponding spins are 
admissible, the Riemannian theta network evaluates to: 

/ a N « 



(36) 



1. 



The degenerate theta network can also be explicitly evaluated; assuming without loss of 
generality that a < b < c: 



(37) 





D 















c> a + b 
c = a + b 
c < a + b. 



Since the Riemannian network's spins are admissible, the third inequality must hold for the 
corresponding areas a, b, c. Thus in this case the conjecture gives an exact formula for the 
Riemannian spin network. 

The next simplest case is the '4j symbol': the spin network with two vertices joined by 
four edges, labelled by positive integers a, b, c and d. Without loss of generality let us assume 
a < b < c < d. As noted in a previous paper in this scries Q, the Riemannian 4j symbol 
counts the dimension of a space of SU(2) intcrtwiners. Using this it follows that: 

R 

'O b + c<d-a 

^{a + b + c — d) d — a<b + c<d + a 
a d + a < b + c. 



(38) 




The corresponding degenerate spin network evaluates to: 



(39) 



so the conjecture is again exact. 





1^ 












d) 



c < d — a 

- a<b + c < d 

- a < b + c, 
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An interesting check on our hypotheses is the tetrahedral spin network. This has four 
vertices and \E\ = ||^|, so the hypotheses of Conjecture 2 do not apply: we expect the 
stationary phase contribution to the Riemannian tetrahedral network to be comparable to 
the degenerate contribution. The degenerate tetrahedral network evaluates to: 

D 



(40) 




247rF(a,5, c, d,e,/)' 



where V{a, b, c, d, e, /) is defined as the volume of the tetrahedron dual to the tetrahedral 
network. Each triangle in this dual tetrahedron corresponds to a vertex of the tetrahedral 
network, and the three sides of the triangle have lengths equal to the labels on the three 
edges incident to the network vertex: 



(41) 



V{a, h, c, d, e, /) = the volume of 




On the other hand, the Riemannian tetrahedral network evaluates to the square of the SU(2) 
tetrahedral network, the basic building-block of the Ponzano-Regge model. Thanks to the 
calculation of Ponzano and Regge later made rigorous by Roberts ||^, this means that: 

R 



(42) 




cos2(5+f) 



'^V2'2'2'2'2'2'' 



_ l + cos2(6'+f) 
oA-jrVlSi k £ i ^ L\' 

zj^ii '^V2'2'2'2'2'2'' 

Here we are dealing with a dual tetrahedron whose edge lengths are ^ati, where Uki ranges 
over a, b, c, d, e, f as 1 < k < I < A. This tetrahedron has Regge action 

l<fc<i<4 

where 9m are the corresponding dihedral angles. Thus it appears that the asymptotics of 
the Riemannian tetrahedral network are a sum of two parts: a part equal to 8 times the 
degenerate tetrahedral network, and an oscillatory part coming from the stationary phase 
points. 

When the edge lengths of the dual tetrahedron are such that it cannot exist in Euclidean 
space, the degenerate tetrahedral network evaluates to zero, and the Riemannian network 
obeys different asymptotics in which its evaluation exponentially decays with increasing spin. 
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6. LORENTZIAN SpIN NETWORKS 

Now let us turn to Lorentzian spin networks ^] . Each positive real number determines 
a unitary irreducible representation of the connnected Lorentz group S0o(3, 1) corresponding 
to an eigenspace of the Laplacian on H^; however, we prefer to describe the evaluation using 
an integral formula. Fixing a connected graph with vertex set V and edge set E, and labelling 
each edge e by a positive real number Oe, we obtain a so-called 'Lorentzian spin network', 
whose evaluation is given by: 

(43) /^a)= [ n<(^(^«(e)'^*(e))) n P^- 

As with degenerate spin networks, we have chosen a vertex vi, set V = V —{vi}, and let x^-^ 
be any fixed point in H^. Here is hyperbolic 3-space, i.e., the submanifold of Minkowski 
spacetime given by: 

^ -x^ -y^ - = 1, t> 0} 
with its induced Riemannian metric. We define the Lorentzian kernel by: 

(44) K^icj^) 



sinh ( 



We warn the reader that this convention differs from that of most previous papers |9[ , 

which include a factor of a in the denominator. Including that factor would divide any 
Lorentzian spin network by the product of its edge labels, so for example, it would divide 
the asymptotics of the Lorentzian lOj symbols as defined here by a factor of A^°. 

The same line of argument by which we arrived at our conjecture concerning asymptotics 
of Riemannian spin networks applies to Lorentzian ones. The most important difference is 
that no factor of 2l^l~^ appears, since there are no 'antipodal points' in hyperbolic space. 
We have not investigated criteria for the existence of stationary phase points, but where they 
are present their exponents will be the same as in the Riemannian case, leading us to make: 

Conjecture 3. Given a connected graph with more than two vertices and \E\ > or 
two vertices and \E\ > 2, or one vertex and \E\ > 1, as A —> oo we have 

/^(Aa)^Al^l-3(l^|-i) /^(a) 

OS long as the integral defining (a) converges and the positive numbers labelling edges 
incident to each vertex are admissible. 

Here we say the positive numbers labelling the edges incident to some vertex are 'admissible' 
if each is strictly less than the sum of the rest. 

Again the simplest test of this conjecture is the Lorentzian theta network. Translating 
their result into our notation, a calculation of Barrett and Crane ||l^ shows that for any 
a,b,c> 0, 



where 



l[/(_a + 5 + c) + /(a-6 + c) + /(a + 6-c)-/(a + 6 + c)] 



/(A;)=tanh(|fc). 
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As the conjecture predicts, the asymptotics of this match those of the degenerate theta 
network, which are given by: 

= - [sign(— a + b + c) + sign(a — 6 + c) + sign(a + b — c) — sign(a + b + c)] . 



It is worth noting that while the integral for the Lorentzian theta network converges even 
when we take the absolute value of the integrand, this fails for the degenerate theta network. 
This makes it more challenging to find criteria for convergence of degenerate spin networks, 
since we cannot simply mimic the theory that applies in the Lorentzian case |p^ . 
Barrett and Crane also worked out the Lorentzian 4j symbols, obtaining: 

^ 7[.9(-a + b + c + d) + g{a ~ b + c + d) + g{a + b ~ c + d) + g{a + b + c - d) 



where 



-g(a + b — c — d) — g{a ~ b + c — d) — g{a — b — c + d) — g(a + b + c + d)], 
g{k) = ^coth(|fc). 



From equation (39) one can show: 



D 



- [h{-a + b + c + d) + h{a-b + c + d) + h{a + b-c + d) + h{a + b + c-d) 
-h{a + b - c - d) - h{a - b + c- d) - h{a - b - c + d) - h(a + b + c + d)], 



where 



Thus the conjecture is also confirmed in this case. 

Next, consider the tetrahedral spin network. As in the Riemannian case. Conjecture 3 
does not apply, but we can still predict the asymptotic behaviour of the contribution of the 
fully degenerate point: 




1 



247rT/(A, A, A, A, A, A) 



An 



SP 
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where 'S'P' represents the contribution from stationary phase points, if any. Below is a 
log-log plot comparing this prediction to numerical data. The horizontal axis in this graph 
represents A, while the vertical axis represents the value of the tetrahedral spin network. 




The most interesting test of Conjecture 3 is the lOj symbol. If the conjecture is true, the 
Lorentzian lOj symbol should be asymptotic to the degenerate lOj symbol, and therefore 
asymptotic to 1/16 times the Riemannian lOj symbol. Since the Riemannian lOj symbol is 
positive 1^, this in turn would imply that the Lorentzian and degenerate lOj symbols are 
positive in the X oo limit. 

It is difficult to compute the Lorentzian lOj symbol, but we have numerically checked the 
conjecture in the special case where all the edges are labelled by the same number A. In this 
case the conjecture says that: 




- .1706 A~^ 

where of course the constant is not exact. Below we show a log- log plot comparing this 
prediction to data points obtained by computing the Lorentzian lOj symbols numerically. 
The horizontal axis represents A, while the vertical axis represents the value of the lOj 
symbol. 
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Here we computed the Lorentzian lOj symbols by applying the VEGAS algorithm to evaluate 
the integral above using polar coordinates on H^. The integral was reduced from 12 to 9 
dimensions by exploiting the S0o(3,l) invariancc; the infinite domain was made compact 
by replacing the radial coordinate for each point, Vi, with a new variable ti = ri/{l + Vi); 
and the domain was further reduced by exploiting a 24-fold symmetry present in the regular 
case. The large dimension and oscillatory nature of this integral make these calculations 
extremely computationally intensive. 

The graph below shows the same Lorentzian 10 j symbols multiplied by A^, plotted on 
a linear scale to give a clearer picture of the rate of convergence towards the asymptote. 
The error bars are three times the standard deviation computed by the VEGAS algorithm. 
The coordinate system we used yielded the lowest standard deviations of several we tried, 
but there is no reason to believe that the estimates of the integral are drawn from a normal 
distribution, so this data should be treated with caution. 
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7. Conclusions 

It appears that the asymptotics of both Riemannian and Lorentzian lOj symbols can 
be computed in terms of degenerate spin networks. For the mathematician, this claim 
still requires proof. Indeed, there is even an issue of rigor concerning our argument that 
degenerate spin networks asymptotically describe the contribution of fully degenerate points, 
since we have not proved that the limit 



v£V' 

exists. The ambitious mathematician could also try to prove more general versions of our 
conjectures, in which all the areas Oe approach infinity, but not in fixed proportion to one 
another. 

For the physicist, however, a more pressing question is: what do these results imply for 
the physics of the Barrett-Crane model? On the one hand, it is unsettling that the simple 
asymptotic behavior of the Ponzano-Regge model is not found here. The way degenerate 
geometries govern the asymptotics of the lOj symbols raises the possibility that in the limit 
of large spins, the Barrett-Crane model reduces to a theory of degenerate metrics. However, 
it is important to bear in mind that the physics of the Barrett-Crane model may not be 
controlled by the large-spin behavior of the lOj symbols — though certainly this affects the 
convergence of the partition function Indeed, there are still major open questions about 
the right way to extract physics from spin foam models, and we hope that our work spurs 
on research on this subject. 

Note Added in Proof. After we submitted this article for publication, preprints of two 
forthcoming papers appeared, one by John Barrett and Chris Steele and another by 
Laurent Freidel and David Louapre ||2^, both of which deal with the asymptotics of lOj 
symbols by somewhat different methods than our own, and both of which obtain results in 
agreement with our conjectures. 
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